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Abstract 

A  computational  model  is  presented  in  which  performance  of  a  solid  oxide  fuel  cell  with  functionally  graded  electrodes  can  be  predicted.  The 
model  calculates  operational  cell  voltages  at  varying  electrode  porosity  and  operational  parameters,  accounting  for  losses  from  mass  transport 
through  the  porous  electrodes,  ohmic  losses  from  current  flow  through  the  electrodes  and  electrolyte,  and  activation  polarization.  Specifically  the 
physical  phenomena  that  occur  when  the  electrode  is  designed  with  a  change  in  microstructure  along  its  thickness  are  studied.  Cell  polarizations  are 
investigated  to  find  arrangements  for  which  the  minimal  polarization  occurs.  Both  diluted  hydrogen  fuel  and  partially  reformed  methane  streams 
with  internal  reforming  are  investigated.  It  is  concluded  that  performance  benefits  are  seen  when  the  electrodes  are  given  an  increase  in  porosity  near 
the  electrolyte  interface  for  certain  fuels  and  with  satisfactory  material  properties.  Enhanced  reformation  is  observed  in  high  tortuosity  structures 
due  to  increased  gas  residence  time  within  the  electrode. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  fuel  flexibility  of  solid  oxide  fuel  cell  (SOFC)  sys¬ 
tems  has  lead  to  their  development  as  next  generation  power 
sources.  The  direct  use  of  liquid  fuels  in  a  single  cell  has  been 
demonstrated  [  1  ] .  The  inherent  high  volumetric  energy  density 
of  hydrocarbon  and  particularly  liquid  fuels  leads  to  theoretical 
increases  in  total  system  power  density.  They  are  therefore  well 
suited  to  a  volume  limited  design  constraint  such  as  air  indepen¬ 
dent  undersea  vehicles  and  other  applications.  Consideration  of 
the  effect  that  this  has  on  the  performance  and  design  of  SOFC 
systems  is  important. 

Functional  grading  of  SOFC  components,  including  the  elec¬ 
trodes  and  electrolyte,  has  been  applied  in  many  successful 
ways.  The  purpose  of  this  can  be  to  increase  electrochemi¬ 
cal  activity  and  to  optimize  the  cell  performance  [1—4],  sepa¬ 
rate  chemically  reactive  components  [3,5],  create  an  interface 
between  two  layers  incompatible  because  of  thermal  expan¬ 
sion  [6],  and  to  increase  mechanical  stability  [7].  Model¬ 
ing  of  the  graded  electrode  manufacturing  process  has  also 
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been  carried  out  to  gain  physical  insight  into  the  graded 
electrode  creation  process  [8,9].  It  is  interesting  to  note  that 
there  have  been  limited  studies  into  modeling  the  mass  trans¬ 
fer  effects  of  electrode  microstructural  changes  during  cell 
performance.  This  paper  addresses  the  mass  transfer  effects 
present  in  the  porous  SOFC  electrodes  with  graded  microstruc¬ 
ture. 

This  study  addresses  modeling  the  performance  of  a  single 
solid  oxide  fuel  cell  operating  on  partially  reformed  methane 
or  diluted  hydrogen  fuel.  It  is  composed  of  an  anode,  elec¬ 
trolyte,  and  cathode  assembly.  The  model  is  able  to  calculate 
voltage  loss  and  species  concentration,  mass  flux,  and  chem¬ 
ical  reaction  rate  profiles.  This  is  achieved  by  including  the 
electrochemical  reactions  at  the  anode/electrolyte  and  the  cath¬ 
ode/electrolyte  interfaces,  internal  steam  reforming  kinetics, 
electric  and  ionic  ohmic  losses,  and  multicomponent  mass  trans¬ 
fer  in  the  porous  electrodes.  This  ID  model  is  suited  for  studying 
the  details  of  the  phenomena  that  occur  in  the  fuel  cell  under 
operating  conditions.  Specifically,  the  changes  that  occur  to  cell 
polarization  and  gas  concentration  distribution  will  be  stud¬ 
ied  to  obtain  physical  insight  and  performance  characteristics 
from  microstructural  features  that  are  both  linearly  changing 
with  depth  and  constant  throughout  the  entire  electrode  thick¬ 
ness. 
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Nomenclature 

Bj  Effective  permeability  coefficient  of  species  i 
(m2  s“ 1 ) 

Cj,  ct  Concentration  of  species  ;,  total  concentration 
(mol  m-3) 

d  thickness  (p,m) 

Dk  Gas  diffusion  constant  (m2  s_  1 ) 

E,  £°  Fuel  cell  electrochemical  potential,  theoretical 
maximum  fuel  cell  electrochemical  potential  (V) 
F  Faraday’s  constant,  96,485  (C) 

i  Current  density  (A  cm-2) 

i0  Exchange  current  density  (A  cm-2) 

K  Knudsen  number  (=)J2r) 

L  Thickness  of  layer  (|xm) 

Nj  Molar  flux  rate  of  species  ;  (mol  m-2  s) 
p  Pressure  (atm.) 

Pi  Partial  pressure  of  gas  (atm.) 

R  Universal  gas  constant,  8.3 1 4  (J  K-moU 1 ) 

(r)  Mean  pore  diameter  (p,m ) 

(r2)  Statistical  distribution  parameter  of  pore  size 
around  mean  (|xm2) 

T  Temperature  (K) 

y;  Concentration  percentage  of  component  i 

Greek  letters 

e  Porosity  of  material 

p  Gas  mixture  viscosity  (m2  s-1) 

Pact  Activation  polarization  (V) 

Pconc  Concentration  polarization  (V) 

Pohm  Ohmic  polarization  (V) 

X  Mean  free  path  length  (m) 

p  Electrical  resistivity  (f2*m) 

r  Tortuosity  of  porous  material 

i j/  Ratio  of  porosity  to  tortuosity  (e/r) 

Subscripts 
A  Anode 

C  Cathode 

E  Electrolyte 

ij  Component  identifier 

Superscripts 
d  Diffusion 

k  Knudsen 

m  Molecular 

op  Operational 

sup  Supply 

TPB  Triple  phase  boundary 

p  Permeation 


2.  Technical  approach 

A  schematic  of  this  SOFC  model  is  shown  in  Fig.  1 .  The  main 
goal  of  the  model  is  to  calculate  the  voltage  of  a  cell  at  a  pre¬ 
scribed  current  draw  as  measured  from  the  surface  of  the  anode  to 


Fig.  1 .  Schematic  of  the  computational  model  domain  and  geometry  of  a  SOFC. 


the  surface  of  the  cathode  for  various  electrode  microstructures. 
AID  model  for  predicting  solid  oxide  fuel  cell  performance  is 
presented  which  will  allow  modeling  of  most  relevant  phenom¬ 
ena. 

The  expression  for  the  operating  voltage  of  the  cell  is  calcu¬ 
lated  according  to: 


Fop  —  F  paC{  7)ohm  Pconc 


0) 


In  Eq.  (1),  Eop  is  the  operational  fuel  cell  voltage,  E  is  the 
ideal  Nernst  voltage  calculated  from  thermodynamic  principles, 
^ohm  is  the  overpotential  caused  by  electronic  and  ionic  ohmic 
losses  from  current  flowing  through  the  fuel  cell  materials,  and 
Pconc  is  the  voltage  loss  caused  by  reduced  fuel  and  oxidant  con¬ 
centrations.  These  are  all  highly  dependent  on  current  density. 

The  Nernst  voltage  can  be  calculated  for  a  H2/O2  fuel  cell 
by: 

n  RT  PU-iPO?0'5 

E  =  E°  + - In  1  2  2  (2) 

2  F  pu2o 

The  thermodynamically  maximum  voltage  obtainable  for  a 
given  temperature  (£°)  is  calculated  according  to  Eq.  (3),  which 
utilizes  the  Gibbs  free  energy  difference  (AG°)  between  fuel  cell 
products  and  reactants,  where  n  =  2  is  the  number  of  equivalent 
electrons  per  mole  reacted.  The  electrochemical  oxidation  of  CO 
is  assumed  to  be  negligible  because  it  is  entirely  consumed  in 
the  reforming  reactions  of  the  cell. 


E° 


—  A  G° 


nF 


(3) 


The  activation  overpotential  can  be  accurately  modeled  by 
the  Butler- Volmer  equation  for  a  fuel  cell  [10]: 


i  —  i0 


2Fpact\ 

RT  J 


—  exp 


— (1  —  >8) 


2Fpact\  3 


(4) 


This  equation  includes  the  exchange  current  density  i0  found 
through  experiments  and  the  transfer  coefficient  [J>  that  is  deter¬ 
mined  through  correlation  with  well  controlled  experiments. 
This  study  assumes  that  functional  grading  of  electrodes,  which 
occurs  in  the  bulk  material,  does  not  change  the  amount  of  elec- 
trochemically  active  sites  at  the  triple  phase  boundary  (TPB). 

The  ohmic  losses  through  the  cell  components,  p0bm,  are 
found  by  using  Ohm’s  law.  The  voltage  drop  is  proportional 
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to  the  current  by  /;0hm.  v  =  iR  where  the  resistance,  R,  is  the  prod¬ 
uct  of  the  material’s  effective  resistivity,  p,  and  the  thickness  of 
the  layer,  L.  For  porous  materials  the  local  effective  resistivity 
is: 


Total  component  flux  can  be  described  by  the  equation: 


P  = 


Pbulk 
1  —  £ 


(5) 


The  total  effective  resistivity  of  an  electrode  is  calculated  by 
finding  the  average  porosity  value  for  a  structure  and  applying 
Eq.  (5).  The  total  ohmic  loss  is  then  the  sum  of  the  three  com¬ 
ponent  losses  from  the  anode,  cathode,  and  electrolyte. 

The  concentration  polarization  occurs  due  to  reduced  reactant 
concentration  at  the  electrochemically  active  regions.  Many  sit¬ 
uations  can  lead  to  significant  resistance  to  mass  transfer  through 
the  porous  electrodes.  This  causes  a  concentration  gradient  of 
products  and  reactants  across  the  electrode  that  leads  to  lower 
concentrations  of  reactants  at  the  triple  phase  boundary  where 
reactants  are  used.  A  dense  layer  of  electrode  is  needed  at  the 
TPB.  This  is  critical  to  increase  the  length  of  the  boundary 
between  the  electrode,  electrolyte,  and  the  gas  phase,  increase 
the  electrochemical  activity  of  the  cell,  and  therefore  enhance 
performance.  It  is  an  assumption  of  the  model  that  the  thickness 
of  the  layer  needed  to  provide  the  enhanced  electrochemical 
activity  was  much  thinner  than  the  anode  thickness  used  for  sup¬ 
port  and  therefore  has  a  negligible  effect  on  the  mass  transfer 
characteristics. 

The  concentration  polarization,  f/cont,  is  the  difference  in 
Nernst  potential  between  the  reactant  supply  channel,  £sup,  and 
the  triple  phase  boundary,  £TPB  [11], 


d  dcx 

Nf  =  -yiBi  —  ,  i=l,...,n 

dx 


(9) 


(r2)Pp  . 


Bi  =  KiV{r)<p  +  (DKiV{r)(  1  -  <f>)  +  =  1,  n  (10) 

8  >1 


<p~l  =  1  +  (  — 
\2(r) 


-l 


(11) 


In  the  above  equations  K,  is  the  Knudsen  number  of  compo¬ 
nent  i,  cu  is  a  physical  factor  that  depends  on  the  Knudsen  wall 
slip,  and  v,  is  the  square  root  of  the  relative  molecular  weights. 
Properties  and  values  for  gasses  used  in  this  study  are  available 
in  literature  [14]. 

The  model  takes  into  account  the  chemical  kinetics  of  both 
the  steam  reforming  and  the  shift  reaction  described  in  Eqs.  (12) 
and  (13). 


CH4  +  H20  <3*  3H2  +  CO 
CO  +  H20  H2  +  C02 


(12) 

(13) 


These  rates  are  determined  experimentally  for  the  general 
temperature  and  pressure  found  within  the  anode  being  mod¬ 
eled  [12,13].  In  order  to  have  a  simple  representation,  they  are 
described  by  a  quasi-homogeneous  expression.  The  rates  of  reac¬ 
tion  of  the  steam  reforming  and  the  shift  reactions  are  described 
by  Eqs.  (14)  and  (15),  respectively. 


>7conc  =  £SUP  -  E 


■TPB 


(6)  Rt  =  kr+ pi P3  -  kr  P2{P4)3 


(14) 


Modeling  mass  transport  in  the  electrodes  is  essential  to 
obtaining  the  necessary  concentrations  for  calculating  the  con¬ 
centration  polarization.  The  mean  transport  pore  model  (MTPM) 
[12,13]  is  used  to  approximate  the  geometry  of  a  porous  media 
as  continuous  round  pores  from  the  electrode  surface  to  the  inter¬ 
face  with  the  electrolyte.  The  relevant  details  are  given  below. 
Within  the  framework  of  the  MTPM,  the  electrode  geometry  is 
described  by  three  parameters  P,  ( r ),  and  (r2). 

The  parameter  \P  can  be  thought  of  as  a  measure  of  the  “effec¬ 
tive”  porosity  of  the  material.  It  adjusts  the  real  porosity  for  the 
extra  resistance  that  the  tortuous  path  imposes  to  the  flow  of 
gasses.  Knudsen  diffusion  must  be  taken  into  account  in  the 
calculations  because  the  pore  sizes  are  on  the  same  order  of 
magnitude  as  the  mean  free  path  length  of  the  gasses.  The  total 
mass  transport  is  given  by: 


Ni  =  N?  +  Nf 


(7) 


where  Nf  is  the  total  component  diffusion  flux  and  Nf  is  the 
total  component  permeation  flux.  This  equation  also  takes  into 
account  Knudsen  diffusion: 


Nf  A  VjNf  -  yiNj  _  dyi 

D\  Dfj  °T  d.x 

7  =  1  1 


(8) 


Rs  —  ks+p2p3  -  ks  P4P5  (15) 

Eqs.  ( 10)— (14)  are  expressions  for  the  molar  rates  of  forma¬ 
tion  of  the  different  species  based  on  the  stoichiometries  of  the 
reactions. 


R<:m4  —  -Rt 

(16) 

Rco  =  R\  —  R& 

(17) 

Rh2o  —  -Rt~  Rs 

(18) 

Rh2  =  37?r  +  Rs 

(19) 

Rco2  —  Rs 

(20) 

For  boundary  conditions,  the  composition  of  the  fuel  gas  is 
known  at  the  interface  between  the  channel  and  the  electrode 
interface.  At  the  interface  between  the  channel  and  electrodes, 
the  molar  flux  rates  of  all  gasses  are  known.  The  resulting  system 
of  ordinary  differential  equations  is  solved  using  an  iterative 
Newton-Raphson  technique.  The  validation  of  the  mass  transfer 
section  of  the  model  is  presented  in  [15].  It  is  also  shown  that  the 
variation  of  P  is  at  least  10  times  more  effective  at  influencing 
the  mass  transfer  in  porous  electrodes  than  other  microstructure 
parameters. 
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3.  Results  and  discussion 

3.1.  Effect  of  the  structural  parameter 

Fig.  2  shows  a  graph  of  numerical  results  of  hydrogen  con¬ 
centration  at  the  triple  phase  boundary  (TPB)  of  a  fuel  cell  anode 
for  different  values  of  the  parameter  constant  across  the  anode 
thickness  and  two  different  fuel  supplies.  The  parameter  has 
a  physical  interpretation  as  the  ratio  of  the  anode  porosity  to  tor¬ 
tuosity.  The  interpretation  of  the  parameter  is  important  through 
much  of  the  paper.  A  structure  with  large  E  can  be  interpreted 
as  either  a  high  porosity  or  low  tortuosity  material.  Both  of  these 
lead  to  a  decrease  in  resistance  to  mass  transfer.  On  the  other 
hand,  a  low  E  leads  to  an  increased  resistance  to  mass  transfer 
because  the  structure  has  low  porosity  or  high  tortuosity. 

The  current  density  is  set  at  500mA cm-2.  Similar  trends 
were  observed  at  other  current  densities.  The  fuel  for  the  dotted 
line  is  20%  H2,  3%  H2O,  and  77%  N2.  For  the  solid  line,  the  fuel 
is  20%  H2  with  25%  CH4, 45%  H20, 7%  CO,  and  3%  C02.  This 
figure  is  an  example  of  the  effect  electrode  microstructure  has 
upon  mass  transfer  within  the  electrodes.  It  also  illustrates  how 
the  incorporation  of  some  amount  of  internal  reforming  vastly 
changes  the  gas  distribution  in  the  anode.  The  concentrations 
are  measured  at  the  TPB  because  this  is  the  concentration  that 
changes  the  operational  voltage  of  a  SOFC.  Hydrogen  at  the  TPB 
falls  significantly  with  an  increase  in  when  internal  reforming 
is  present.  When  a  diluted  hydrogen  stream  is  fed  to  the  fuel  cell, 
the  partial  pressure  of  H2  increases  with  an  increase  in  the  value 
of  the  E  parameter.  This  can  be  explained  by  the  source  of  H2 
that  internal  reforming  introduces  into  the  system.  The  longer 
residence  time  of  CH4  in  the  anode  microstructure,  caused  by 
the  lower  E  value  (more  tortuous  channels),  leads  to  an  increase 
in  H2  production  within  the  electrode  volume.  These  results 
strongly  suggest  that  microstructural  grading  of  the  electrode 
can  improve  fuel  cell  performance. 

3.2.  Polarization  effects 

Section  3.1  has  shown  that  electrode  microstructure  for  dif¬ 
ferent  fuel  compositions  has  substantial  impact  on  cell  perfor- 


H* 

Fig.  2.  Numerical  predictions  of  *1*  vs.  partial  pressure  of  H2  at  the  TPB  for 
a  current  density  of  500  m.A  cm  2.  Results  show  that  this  value  changes  very 
differently  with  'f'  depending  on  fuel  composition. 


Current  Density  (mA/cm2) 

Fig.  3 .  Total  cell  voltage  for  a  20%  H2/3%H2  0/77 %  N2fuel  stream  with  different 
anode  microstructures.  Maximum  power  density  is  8.7%  higher  for  case  {4}  than 
for  case  {1}. 

mance.  This  is  further  shown  through  polarization  curves  in 
Figs.  3  and  4,  where  graphs  of  voltage  versus  current  density  are 
given  for  a  variety  of  constant  and  graded  microstructures.  Four 
cases  are  presented  and  used  henceforth.  Values  of  E  from  0.05 
to  0.25  are  modeled,  approximating  minimum  and  maximum 
physically  realistic  values  from  literature  [8,12],  respectively. 
The  cases  are: 

•  {1}  and  {2}  are  when  the  E  parameter  is  kept  constant  at  a 
representative  maximum  and  minimum  value  (0.05  and  0.25) 
across  the  anode  thickness,  respectively. 

•  {3}  is  when  the  parameter  across  the  anode  thickness  is 
varied  linearly  from  0.25  at  the  fuel  inlet  channel  to  0.05  at 
the  electrolyte  interface. 

•  {4}  is  when  the  parameter  across  the  anode  thickness  is 
varied  linearly  from  0.05  at  the  fuel  inlet  channel  to  0.25  at 
the  electrolyte  interface. 

The  numerical  model  showed  that  the  resulting  cases  in 
between  these  extreme  values  of  predictably  lie  within  the 
results  given  by  these  values.  Linear  grading  of  SOFC  materials 
has  been  obtained  using  advanced  ceramic  techniques  [7],  An 
electrode  can  also  be  composed  of  layers  of  constant  (psi).  In 
this  study,  the  layers  of  constant  will  be  represented  as  a  linear 
variation  along  the  anode  thickness. 


Current  Density  (mA/cm2) 

Fig.  4.  Total  cell  voltage  for  a  20%  H2/25%  CH4/45%  H20/7%  CO/3%  C02 
fuel  stream  with  various  alternative  anode  microstructures.  Maximum  power 
density  is  19%  higher  for  case  {2}  than  for  case  {1}. 
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Table  1 

List  of  parameters  used  in  this  study 


Parameter 

Value 

Temperature  (T,  K) 

1273 

Pressure  ( P ,  atm) 

1.00 

Exchange  current  density  (i0,  A/cm2) 

0.05 

Transfer  coefficient  (f3) 

0.5 

Anode  thickness  (La,  mm) 

1.0 

Anode  bulk  resistivity  (p,  £2/m) 

0.001 

Anode  tortuosity  (ta) 

2.00 

Anode  pore  radius  (r,  |xm) 

1.05 

Cathode  thickness  (Lc,  |xm) 

20 

Cathode  bulk  resistivity  (p,  £2/m) 

3.5 

Cathode  tortuosity  (rc) 

2.00 

Cathode  pore  radius  ((/*),  [xm) 

0.75 

Electrolyte  thickness  (Le,  |xm) 

10 

Electrolyte  resistivity  (p^,  Q/m) 

0.20 

y  Profile 
{1)0.25  ->0.25 
{2}  0.05  ->  0.05 
13}  025"T0.05 
14}  005^0.25 


Fig.  5.  Hydrogen  concentration  profiles  through  the  anode  thickness  for 
microstructure  cases  { 1 } — {4}  with  diluted  hydrogen  fuel.  The  TPB  is  located 
at  1  mm  depth. 


In  Fig.  3  the  fuel  stream  provided  to  the  cell  is  20%  H2,  3% 
H2O,  and  77%  No.  Fig.  4  shows  an  internal  reforming  SOFC 
with  a  20%  H2,  25%  CH4,  45%  H20,  7%  CO,  and  3%  C02  fuel 
stream.  Other  parameters  considered  in  this  study  are  shown  in 
Table  1.  With  these  cases,  the  effect  of  grading  the  microstruc¬ 
ture  of  a  SOFC  anode  upon  the  actual  fuel  cell  performance  is 
predicted. 

In  Fig.  3,  minimal  polarization  is  observed  to  occur  for  case 
{4},  which  has  a  graded  microstructure.  Only  a  slightly  better 
performance  is  seen  for  this  case  than  for  no  microstructure  grad¬ 
ing  (case  {2}).  The  ohmic  voltage  losses  are  also  shown  in  this 
figure.  Both  electrode  grading  methods  produced  a  microstruc¬ 
ture  with  the  same  average  resistance,  as  expected.  The  lowest 
ohmic  voltage  loss  is  seen  for  case  {2}  because  low  porosity 
material  provides  more  material  for  electronic  conduction.  The 
worst  ohmic  voltage  loss  is  seen  for  the  high  d*  electrode  {1} 
due  to  its  high  resistance.  And  the  graded  electrodes  have  a  per¬ 
formance  in  between  these  two  constant  microstructure  cases, 
and  do  not  reflect  the  same  trend  shown  in  the  polarization  curve 
since  concentration  and  activation  polarization  also  plays  a  role 
in  cell  performance.  We  observed  the  lowest  concentration  polar¬ 
ization  at  maximum  power  for  cases  {1}  and  {4}  (results  not 
shown).  This  can  be  attributed  the  high  species  counter-diffusion 
flux  area  at  the  TPB  where  H20  is  produced  and  H2  is  actively 
consumed,  producing  large  concentration  gradients.  The  open 
pore  structure  at  the  TPB  given  by  the  graded  electrode  allows 
more  effective  gas  transport  while  still  keeping  ohmic  losses 
low.  It  is  important  to  note  that  current  cell  designs  have  outer 
layers  that  are  relatively  porous  with  a  thin  low  porosity  layer 
near  the  TPB.  The  philosophy  behind  this  is  to  not  limit  mass 
transfer  with  a  thick  layer  of  low  porosity  material  near  the  chan¬ 
nel  and  to  increase  the  active  TPB  area.  This  modeling  suggests 
that  alternative  electrode  structural  designs  are  worth  consider¬ 
ing.  It  should  also  be  noted  that  most  functional  grading  methods 
have  a  performance  that  is  minimally  lower  than  a  homogeneous 
electrode  for  the  H2  fed  fuel  cell. 

Fig.  4  displays  the  additional  effect  of  internal  reforming  in 
the  anode.  The  Nernst  voltage  was  calculated  with  the  values 
of  H2  and  H20  concentrations  at  the  TPB.  The  concentration 
polarization  does  not  have  a  meaningful  reference  in  this  case 


because,  due  to  internal  reforming,  the  channel  concentration  of 
H2  does  not  represent  a  significant  value  for  comparison.  The 
result  is  markedly  different  from  Fig.  3.  The  best  performance  is 
seen  for  case  {2},  which  represents  a  homogeneous  low  poros¬ 
ity/high  tortuosity  structure. 

This  means  that  there  are  limited  benefits  seen  under  these 
operating  conditions  to  justify  the  extra  cost  of  manufactur¬ 
ing  the  graded  structure.  The  internal  reforming  of  CH4  leads 
to  a  large  increase  in  H2  concentration  throughout  the  anode, 
from  10%  at  the  inlet  to  typically  30%  at  the  TPB.  In  fact,  the 
efficiency  of  this  reforming  process  is  a  large  driver  for  the  per¬ 
formance  of  the  cell.  It  is  also  observed  that  one  of  the  main 
reasons  for  the  increase  in  performance  for  case  {2}  is  that  it 
has  the  lowest  ohmic  losses.  We  observed  that  the  key  to  the 
better  performance  is  an  effective  reforming  reaction.  If  lower 
levels  of  reforming  are  seen,  the  system  will  begin  to  resemble 
a  humidified  H2  fuelled  system  and  performance  may  be  gained 
by  electrode  grading.  An  analysis  of  the  fuel  cell  operating  range 
would  be  required  to  make  a  decision  on  the  suitability  of  this 
scheme. 

The  effect  of  changing  the  cathode  porosity  was  also  studied. 
It  was  computed  that  the  same  variation  of  the  microstruc- 
tural  parameter  as  given  for  cases  { 1 }— {4}  above  leads  to 
a  variation  in  performance  of  less  than  1%.  This  is  attributed 
to  the  lack  of  counter-diffusion  that  occurs  in  this  electrode 
as  the  oxygen  ions  are  directly  transported  through  the  elec¬ 
trolyte.  This  observation  also  suggests  the  possible  application 
of  functionally  graded  electrodes  to  anionic  conducting  elec¬ 
trolytes. 

3.3.  Species  diffusion  in  graded  microstructure  electrodes 

The  species  within  the  fuel  cell  are  investigated  at 
500  mA  cm-2  current  density  and  conditions  described  in 
Table  1 .  The  figures  will  provide  insight  into  the  physical  mech¬ 
anisms  that  drive  the  changes  in  performance  that  are  seen  in 
Figs.  3  and  4.  Hydrogen  will  be  the  species  studied  because  of 
its  obvious  importance  to  fuel  cell  performance. 

Fig.  5  is  a  graph  of  the  partial  pressure  of  H2  throughout 
the  depth  of  the  anode  for  the  four  cases  of  grading  when 
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Fig.  6.  (a)  Hydrogen  flux  profiles  through  the  anode  thickness  for  microstructure 
cases  { 1 } — [4}  with  internal  reforming.  The  TPB  is  located  at  1  mm  depth,  (b) 
hydrogen  concentration  profiles  through  the  anode  thickness  for  microstruc¬ 
ture  cases  { 1  }-{4}  with  internal  reforming.  The  TPB  is  located  at  1  mm 
depth. 


the  SOFC  is  fed  with  diluted  H2  fuel.  It  is  seen  that  the  same 
concentration  profile  is  predicted  for  each  case  where  the  value 
at  the  TPB  is  the  same.  This  is  evidence  of  the  importance  of 
the  TPB  upon  the  mass  transfer  in  a  operational  SOFC.  Because 
the  fuel  cell  is  operating  at  steady  state  it  can  only  maintain  a 
linear  concentration  profile  when  there  are  no  internal  sources  of 
species,  as  is  the  case  of  a  diluted  hydrogen  fuel.  This  leads  to  the 
concentration  profile  being  determined  by  the  largest  resistance 
in  the  electrode,  which  is  the  area  at  the  TPB.  A  high  'P  value 
at  the  TPB  enhances  mass  transfer  in  this  area  and  therefore 
increases  the  concentration  of  FI2  at  the  TPB. 

Evidence  for  the  effect  of  counter-diffusion  of  species  can 
be  seen  in  the  H2  species  molar  flux  profiles  that  are  shown  in 
Fig.  6(a).  These  cases  are  all  for  identical  gas  inlet  conditions 
described  in  Table  1  with  internal  reforming.  The  molar  flux  rate 
of  Fh  is  shown  through  the  thickness  of  the  anode.  The  TPB  is  at 
a  depth  of  1  mm.  Hydrogen  gas  is  volumetrically  produced  in  the 
anode  from  methane  and  consumed  at  the  TPB.  It  is  generally 
seen  that  the  flux  profile  is  most  strongly  dependent  upon  the 
>P  value  that  occurs  near  the  TPB.  The  greatest  negative  values 
of  H2  flux  are  seen  for  case  {4}.  For  case  {1},  where  there  is 
also  the  large  pore  structures  near  the  TPB,  the  H2  flux  is  also 
significantly  more  negative.  These  results  once  again  support  the 
conclusion  that  an  open  pore  structure  near  the  electrochemically 
active  layer  will  enhance  the  counter  diffusion  of  species  that 
must  occur  in  this  region. 


This  phenomena  is  also  visible  is  the  concentration  profiles  as 
seen  in  Fig.  6(b),  which  shows  hydrogen  concentration  through 
the  anode  thickness.  This  graph  shows  a  significant  increase 
in  the  H2  partial  pressure  for  the  case  where  the  microstruc¬ 
ture  is  graded  from  low  *P  at  the  fuel  channel  to  high  at  the 
electrolyte  interface.  However,  it  was  shown  in  Fig.  4  that  this 
increase  in  electrochemical  performance  is  negated  by  the  reduc¬ 
tion  in  electrical  conductivity.  Anode  materials  with  increased 
conductivity  could  change  this  balance.  It  can  also  be  seen  in 
Fig.  6(b)  that  there  are  significant  changes  to  the  rate  of  change 
of  hydrogen  concentration  throughout  the  microstructure  thick¬ 
ness.  It  is  clear  that  case  {4}  produces  a  sharp  rise  in  hydrogen 
concentration  that  is  triggered  by  the  residence  time  of  CH4  in 
the  surface  of  the  anode.  It  appears  that  this  is  actually  enhanced 
by  the  increased  mass  transfer  rate  that  occurs  because  of  the 
lowering  of  the  'P  parameter  throughout  the  anode  thickness 
that  allows  the  counter  diffusion  of  reactants  to  occur  near 
the  TPB.  Furthermore,  it  is  noted  that  the  smallest  TPB  H2 
concentration  is  created  by  case  {1},  which  is  also  the  worst 
electrical  conductor.  This  leads  to  its  theoretically  low  perfor¬ 
mance. 

4.  Conclusions 

The  effect  of  electrode  microstructure  grading  on  the  per¬ 
formance  of  a  solid  oxide  fuel  cell  has  been  studied.  This  was 
performed  to  determine  possible  ways  that  the  performance  can 
be  improved  with  a  simple  design.  This  goal  was  accomplished 
by  utilizing  a  numerical  model  of  the  fuel  cell  electrodes,  elec¬ 
trolyte,  and  gas  channel  to  predict  performance  for  alternative 
electrode  geometries. 

It  was  shown  that  when  diluted  H2  is  the  fuel  stream,  grad¬ 
ing  the  electrode  from  less  to  more  porous  towards  the  TPB, 
increases  performance  beyond  having  a  monolith  of  constant 
1 P  of  either  extreme  value.  Also,  it  was  observed  that  when  the 
microstructure  is  graded  for  an  internal  reforming  fuel  cell  there 
are  no  performance  gains  when  available  materials  are  used.  It 
was  observed  that  higher  concentrations  of  H2  were  seen  at  the 
TPB  when  an  optimal  microstructure  was  used  for  an  internal 
reforming  fuel  cell,  but  ohmic  losses  negated  this  advantage.  A 
detailed  investigation  of  non-isothermal  effects  is  necessary  to 
determine  the  interplay  of  heat  and  mass  transfer,  and  its  effect 
on  electrode  grading.  The  use  of  an  electrode  that  is  graded 
from  high  *P  at  the  TPB  to  low  {P  at  the  supply  channel  has  been 
shown  to  be  an  effective  way  to  promote  diffusion  mass  transfer 
when  diluted  hydrogen  in  introduced  as  a  fuel  to  a  SOFC  with 
a  significantly  thick  anode. 
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